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Within the picture of jet quenching induced by multiple parton scattering and gluon 
bremsstrahlung, medium modification of parton fragmentation functions and therefore the sup- 
pression of large transverse momentum hadron spectra are controlled by both the value and the 
space-time profile of the jet transport parameter along the jet propagation path. Experimental 
data on single hadron suppression in high-energy heavy-ion collisions at the RHIC energy are ana- 
lyzed within the higher-twist (HT) approach to the medium modified fragmentation functions and 
the next-to-leading order (NLO) perturbative QCD (pQCD) parton model. Assuming that the jet 
transport parameter q is proportional to the particle number density in both QGP and hadronic 
phase, experimental data on jet quenching in deeply inelastic scattering (DIS) off nuclear targets 
can provide guidance on qh in the hot hadronic matter. One can then study the dependence of 
extracted initial value of jet quenching parameter go at initial time to on the bulk medium evolu- 
tion. Effects of transverse expansion, radial flow, phase transition and non-equilibrium evolution 
are examined. The extracted values are found to vary from <jot~o = 0.54 GeV 2 in the (l+3)d ideal 
hydrodynamic model to 0.96 GeV 2 in a cascade model, with the main differences coming from the 
initial non-equilibrium evolution and the later hadronic evolution. The overall contribution to jet 
quenching from the hadronic phase, about 22-44%, is found to be significant. Therefore, realistic 
description of the early non-equilibrium parton evolution and later hadronic interaction will be crit- 
ical for accurate extraction of the jet transport parameter in the strongly interacting QGP phase in 
high-energy heavy-ion collisions. 

PACS numbers: 12.38. Mh,24.85.+p,25.75.-q 



I. INTRODUCTION 



A large amount of experimental data from the Rela- 
tivistic Heavy Ion Collider (RHIC) PHI] strongly sug- 
gest that a novel form of matter, a strongly coupled 
quark gluon plasma (sQGP), may have been formed in 
the central region of high-energy heavy-ion collisions. 
Among the experimental evidences for the formation of 
sQGP arc the jet quenching phenomena that include the 
strong suppression of hadron spectra with large trans- 
verse momentum in central Au + Au collisions as com- 
pared to p + p collisions [!, and the disappearance of 
back-to-back correlations of large transverse momentum 
hadrons 0]. These jet quenching patterns observed at 
RHIC arc in good agreement with the theoretical pre- 
dictions for jet quenching (l-fil^ or suppression of large 
transverse momentum hadrons as a consequence of par- 
ton energy loss and medium modified parton fragmenta- 
tion functions induced by multiple scattering as partons 
propagate through the dense medium after their initial 
production. The parton energy loss and medium mod- 
ification of the fragmentation functions due to multiple 
parton scattering and induced gluon bremsstrahlung arc 



controlled by the jet transport parameter [14j . 

QR = P j d QT C -^T<lT- (!) 

or the average squared transverse momentum broadening 
per unit length for a jet in color representation R, which 
is also related to the gluon distribution density of the 
medium [3, [HI and therefore characterizes the medium 
property as probed by an energetic jet. Here we consider 
a picture in which the jet parton interacts with medium 
particles or quasi-particlcs with cross section an and p is 
the local particle density of the medium. Study of the jet 
quenching phenomena therefore can provides important 
information on the space-time profile of the jet transport 
parameter and consequently properties of the sQGP in 
heavy-ion collisions. 

Extensive phenomenological studies on the suppres- 
sion of single hadron spectra [HI, Illl|l9j and dihadron 
correlations [2(| [2l| at large transverse momentum have 
been carried out since the first observation of the strong 
jet quenching phenomena. Recent emphasis of the phe- 
nomenological studies has been shifted to systematic 
analyses of the experimental data with different jet 
quenching models and qualitative extraction of the jet 
transport parameter [22h26| . In this paper, we will carry 
out phenomenological analysis of the suppression of sin- 
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FIG. 1: A typical Feymann diagram for quark-gluon re- 
scattering processes in DIS with three possible cuts, cen- 
tral^), left(L), and right(R). 



II. ENERGY LOSS AND MODIFIED 
FRAGMENTATION FUNCTIONS 

Within the generalized factorization of twist-four pro- 
cesses, one can calculate the nuclear modification of frag- 
mentation functions and the energy loss of a quark prop- 
agating through nuclear matter after it is produced via a 
hard process in DIS off a nuclear target [2^, [29[ . Within 
such HT approach, the medium modification to the quark 
fragmentation functions in DIS off a nuclear target is 
caused by multiple scattering between the struck quark 
and the nuclear medium on the quark's propagation path 
and the induced gluon bremsstrahlung as illustrated in 
Fig. ([1]). The medium modified quark fragmentation 
function, 



gle hadron spectra within the higher-twist (HT) formal- 
ism for medium modification of the parton fragmenta- 
tion functions and the next-to-leading order perturbative 
QCD parton model [2(| for initial jet production. We 
will calculate the single hadron suppression within three 
different models for the dynamical evolution of the bulk 
matter: (l+l)d Bjorken model, (l+3)d ideal hydrody- 
namic model and a parton cascade model. We will study 
in particular the effect of collective expansion, transverse 
flow, non-equilibrium and most importantly jet quench- 
ing in the hadronic phase of the bulk matter. Since the 
jet transport parameter is related to the gluon distribu- 
tion density of the bulk medium, one expects it to be 
different in a QGP and hadronic matter. We will use the 
information on jet transport parameter in cold nuclei ex- 
tracted from phenomenological studies of deeply inelastic 
scattering (DIS) off large nuclei [13] and extrapolate to a 
hot hadronic gas assuming that the gluon distribution in 
a nucleon in cold nuclei is the same as in a hot hadronic 
gas. This allows us to focus on the values of the jet trans- 
port parameter in a QGP that can be accommodated by 
the experimental data on single hadron suppression and 
the effect of dynamic evolution of the bulk medium. 



The rest of the paper is organized as follows. In the 
next section, we will review the energy loss and modified 
fragmentation functions in hot nuclear medium within 
the HT approach and their dependence on the space- 
time profile of the jet transport parameter qj^. We then 
describe the next-to-leading order (NLO) pQCD parton 
model for single hadron production in heavy-ion colli- 
sions in Sec. IIIII In Sec . IIV( we will introduce three 
different dynamic evolutions for the bulk medium. The 
numerical calculations and phenomenological studies of 
the experimental data on single hadron suppression and 
extraction of the jet transport parameters within each 
dynamic model are carried out in Sec . [V] Finally, we 
conclude in Sec. I VII with a summary. 
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takes a form very similar to the vacuum bremsstrahlung 
corrections that leads to the Dokshitzer-Gribov-Lipatov- 
Altarclli-Parisi (DGLAP) [3(3] evolution equations in 
pQCD for fragmentation functions, except that the 
medium modified splitting functions, 
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depend on the properties of the medium through the 
twist-four quark-gluon correlations inside the nucleus, 
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where £ = y\ — y%, y, Hi and yi are space-time coordinates 
associated with the quark and gluon fields as illustrated 
in Fig. [TJ Here we include contributions beyond the he- 
licity approximation [2^ | and consider only the twist-four 
matrices that are enhanced by the large nuclear size. The 
relative transverse coordinate £t is the Fourier conjugate 
of the transverse momentum qx in the gluon distribution 
function. The momentum fraction 
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FIG. 2: Feymann diagram for induced gluon radiation in hot 
QGP medium that contributes the quark energy loss. 



is the total (+component) longitudinal momentum trans- 
fer from the target to the propagating quark and radiated 
gluon with longitudinal momentum zq~ and transverse 
momentum if. We impose kinematic constraints xl < 1 
and £t < min[zE, (1 — z)E) in the integration in Eq. ([2]). 
A similar (+) longitudinal momentum transfer 
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is also provided by the initial gluon with transverse mo- 
mentum <jt- In the modified splitting function, 



^le^ + y- (A\^ q (0)^ q (y-)\A) (8) 

Z7T Z 



is the quark distribution function which represents the 
production rate of the initial quark carrying x = 
Q 2 /2p + q~ (the Bjorken variable) fraction of the nucleon 
(+) longitudinal momentum in DIS. The quark-gluon 
matrix element, 
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(9) 



comes from the virtual correction to the induced gluon 
bremsstrahlung process. 

If we define the quark energy loss as the energy car- 
ried by the radiative gluon in the multiple scattering pro- 
cesses, the total energy loss for a propagating quark in a 
deeply inelastic scattering (DIS) off a large nucleus is 
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One can extend the results for medium modified parton 
fragmentation functions in DIS to the case of quark prop- 
agation in a hot medium (either QGP or hot hadronic 
matter) after it is produced through initial hard processes 
before the formation of the QGP [JJ|, [22[ in high-energy 
heavy-ion collisions. In this case, we can assume the 
quark is produced through a hard process such as the 
virtual-photon- nucleus collisions in Fig. [5] and will inter- 
act with partons in the hot medium which is independent 



of the initial quark production process. Therefore, we 
can replace the nucleus state as a product of the initial 
nucleus state and a thermal ensemble of quasi-particle 
states in the hot medium, 
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and consider only the final state interaction between the 
produced quark and the bulk medium as illustrated in 
Fig. [2J Here f(k,y) = \ipk(y)\ 2 is the local phase-space 
density of the quasi-particle distribution in the medium. 
We also neglect multiple particle correlations inside the 
hot medium which can be included through an effective 
gluon distribution density [l5j . Under these assumptions, 
the quark-gluon correlation function in the HT approach 
to multiple scattering will take a factorized form, 
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where c(x,xl) = fq{x + xi,)/fq{x) and q R is the gener- 
alized jet transport parameter, 
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which depends on the fractional (+) longitudinal momen- 
tum transfer xl from medium gluons [3l|. Here <fik(x, qr) 
is the transverse momentum dependent gluon distribu- 
tion function from the quasi-particle in the medium, 
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The two terms in the square brackets of Eq. ([12]) cor- 
respond to two different processes in the multiple scat- 
tering. The first term involves (+) longitudinal momen- 
tum transfer xl between the propagating parton and the 
medium due to final gluon production. It contains what 
is normally defined as pure elastic energy loss [32j . The 
second term that is proportional to the normal (or spe- 
cial) transport parameter q R (E,y) = q R (E,XL = 0,y) 
corresponds to pure radiative processes. In this paper, 
we assume x 3> xl , xt and only focus on the contribu- 
tion of radiative energy loss. Therefore, c{x,xl) ~ 1, 
q R {E,x L ,y) w q R (E, 0, y) = q R (E,y). Given the space- 
time profile of the jet transport parameter q R (E, y), one 
should be able to calculate the modified fragmentation 
function according to Eq. The corresponding quark 
energy loss in Eq. ([TO]) can be expressed as 
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in terms of the jet transport parameter. The jet trans- 
port parameter for a gluon is 9/4 times of a quark and 
therefore the radiative energy loss of a gluon jet is also 
9/4 times larger than that of a quark jet. In our following 
calculations of the hadron spectra in heavy-ion collisions, 
we will use the same formalism for both quark and gluon 
modified fragmentation functions but with jet transport 
parameters that differ by a factor of 9/4 for quark and 
gluon jets. 



III. SINGLE HADRON SPECTRA WITHIN NLO 
PQCD PARTON MODEL 

In this paper, we employ the NLO pQCD parton model 
for the initial jet production spectra which has been 
shown to work well for large pr hadron production in 
high energy nucleon-nucleon collisions [331 ]. In a factor- 
ized form, the inclusive particle production cross section 
in p + p collisions can be calculated as a convolution of 
parton distribution functions inside the proton, elemen- 
tary parton-parton scattering cross sections and parton 
fragmentation functions, 
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where da/dt(ab — > cd) are elementary parton scatter- 
ing cross sections at leading order (LO) a 2 . The next- 
to-leading order (NLO) contributions including 2^3 
tree level contributions and 1-loop virtual corrections 
to 2 — > 2 tree processes. We refer to Ref. [20| for 
a detailed description. The proton parton distribution 
functions (PDF) L/ P (x a ,fi 2 ) are given by the CTEQ6M 
parametrization |34| where x a is the fractional momen- 
tum carried by the beam partons. The fragmentation 
function (FF) D^ c (z c , /i 2 ) of parton c into hadron h, 
with z c the momentum fraction of a parton jet carried 
by a produced hadron are given by the updated AKK 
parametrization (35[, which is recently improved with 
new input from the RHIC data, especially for K ± , 
p/p, Kg and A/A particles. 

We use a NLO Monte Carlo based program [20| to cal- 
culate the single hadron spectra in our study. In this 
NLO program, the factorization scale and the rcnormal- 
ization scale are chosen to be the same (denoted as fi) and 
are all proportional to the transverse momentum of the 
final hadron pt- As pointed out in Ref. [2lJ, the calcu- 
lated single inclusive pion spectra within the NLO pQCD 
parton model for p + p collisions agree well with the ex- 
perimental data at the RHIC energy using the scale in 
the range fj, = 0.9 ~ l.bpx- We find that the calculated 
7T spectra in p + p collisions with the scale /i = \.2pr 
can fit RHIC data very well and will use the same scale 
in A + A collisions in our following calculations. 



We further assume that the factorized form for large 
transverse momentum single hadron spectra in NLO 
pQCD parton model can be applied to heavy-ion colli- 
sions, 
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where (b) is the number of binary nucleon-nucleon 
collisions at a fixed impact b of A + B nuclear collisions, 
(D^(z h ,Q 2 ,E 7 b)} is the medium modified parton frag- 
mentation function averaged over the initial production 
position and jet propagation direction, and f a /A(x a , /J 2 ) 
is the effective parton distributions per nucleon inside 
a nucleus. For large transverse momentum hadron pro- 
duction, the effect of nuclear modification of the parton 
distributions is small. We will neglect such nuclear effect 
in our study here. 

According to the definition of the jet transport param- 
eter in Eq. (|13p . it should be proportional to the local 
particle density at fj = f + (t — r^jn in either QGP or 
hadronic phase of the evolving bulk medium along the 
path of a propagating jet which is a straight line in the 
cikonal approximation. Here tq is the formation time of 
the medium and the direction of jet propagation is de- 
fined by its azimuthal angle </> in the transverse plane. 
The quark-gluon correlation in Eq. (fT2"j) that enters the 
modified fragmentation function in Eq. ([5]) will involve an 
integration over the duration of the jet propagation and 
depend on the initial production position f and propaga- 
tion direction ft. In heavy-ion collisions at a fixed impact 
parameter b, the initial jet production cross section is 
proportional to the overlap function of two colliding nu- 
clei J d 2 rtA(r)tB(\b — r\). One therefore has to average 
over the initial jet production position and propagation 
direction to obtain the effective medium modified parton 
fragmentation function, 
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With a given space-time profile of the jet transport 
parameter, one can therefore use the above effective 
medium modified fragmentation functions to calculate 
the large transverse momentum hadron spectra and the 
the suppression factor (or nuclear modification factor), 
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for characterization of the effect of jet quenching on 
hadron spectra in heavy-ion collisions [361 ] . Phenomeno- 
logical study of the experimental data on the above 
hadron suppression factor will in turn provide us con- 
straints on the space-time profile of the jet transport pa- 
rameter <jr(E, y). 

The computation of the effective medium-modified 
jet fragmentation function in Eq. (|18[) involves a 6- 
dimensional integration which is in addition to the Monte 
Carlo integration in the NLO pQCD calculation of the 
single hadron spectra [2(| ■ Numerically, we will compute 
and tabulate the medium modified fragmentation func- 
tions as functions of the initial jet energy E, fractional 
momentum z and the factorization scale Q 2 for fixed 
value of impact parameter b. We will then use numerical 
interpolation of the table for the medium modified frag- 
mentation functions in the NLO pQCD calculation of the 
single hadron spectra. 



IV. BULK MATTER EVOLUTION 



Recent phenomenological studies [15j, |22| - |25| have fo- 
cused on the initial values of the jet transport parame- 
ters as constrained by the experimental data. We study 
in this paper the uncertainty of the extracted initial jet 
transport parameter due to different models of the dy- 
namical evolution of the bulk medium and in particu- 
lar the effect of transverse expansion, flow effect, non- 
equilibrium and phase transition by considering three 
different models for the bulk evolution. Assuming fast 
thermalization of the partonic matter in the initial stage 
of high-energy heavy-ion collisions, we will consider first 
(l+l)d Bjorken ^ and (l+3)d [M [I| ideal hydro- 
dynamic evolution of the bulk matter for the calcu- 
lation of medium modified parton fragmentation func- 
tions. Inclusion of shear and bulk viscosity in the vis- 
cous hydrodynamics can affect significantly the elliptic 
flow or azimuthal asymmetry of the final hadron spectra 
at freeze-out. Their effect on the space-time profile of the 
bulk medium such as entropy density, however, is much 
smaller (on the order of a few percent) pol - l43l | which we 
will neglect here in the hydrodynamic evolution of the 
bulk matter for the study of jet quenching. To consider 
non-equilibrium evolution of the bulk matter and its ef- 
fect on jet quenching, we will also calculate medium mod- 
ified hadron spectra with the space-time profile of parton 
density as given by a parton cascade model [3 145| ■ 



A. (l + l)d Bjorken Expansion 

In a (l+l)d Bjorken model [37j for the bulk evolution, 
the system is assumed to undergo a longitudinally boost- 
invariant expansion with the (l+l)d ideal hydrodynamic 
equation, 



de 



P 



(20) 



where e is the energy density, P is the pressure and 
r = \/t 2 — z 2 is invariant time. For a massless ideal gas 
equation of state (EOS), the solution of the above ideal 
hydrodynamic equation leads to a time evolution of the 
entropy density, 



(21) 



with s the entropy density at the initial time r . Since 
the jet transport parameter qji is directly proportional to 
gluon distribution density and therefore proportional to 
the particle number or entropy density, we assume that 
it will have a similar time dependence due to longitudi- 
nal expansion. In the (l+l)d Bjorken model, we further 
neglects the transverse expansion and phase transition. 
Therefore the transverse profile of jet transport param- 
eter will be given by the initial transverse profile of the 
parton density which we assume to be proportional to 
the transverse density of participant nucleons according 
to the wounded nucleon model of nucleus-nucleus colli- 
sion 14611. 
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where n par t (6, r) is the transverse density of participant 
nucleons in nucleus-nucleus collisions with impact param- 
eter b, 



''part 



and ctnn is the nucleon-nucleon total inelastic cross sec- 
tion, which is set to be <tnn = 42 mb at the RHIC energy. 
Here t A {f) is the nuclear thickness function, 
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and p A (f,z) is the single nucleon density in the nu- 
cleus normalized to J d 2 rt A (f) = A. The total number 
of participant nucleons in nucleus-nucleus collisions at 
impact-parameter b is then N petrt {b) = / d 2 rn par t(6, r). 
In the space-time profile of the jet transport parameter 
in Eq. (|22[) . qo is defined as the jet transport parame- 
ter at the center of the bulk medium in the most central 
nucleus-nucleus collisions (b = 0). 

We will consider two different nuclear distributions for 
p A (r, z) in the (l+l)d Bjorken model. One is the Woods- 
Saxon distribution: 



PA{r,z) 
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where d = 0.662 fm, R A = 1.26^4 1 / 3 fm, and no is the 
normalization factor given by J d 2 rdz p A (r, z) = A pTl ]. 
We also consider the simple hard-sphere distribution, 
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with Ra = 1.12A 1 / 3 for comparison. 

In the present work we focus on jet quenching in the 
10% most central collisions. By comparing the averaged 
number of participating nuclcons, 



(N pait } 



J dbbN paxt (b) 
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which is given by experiments for the most 10% cen- 
tral collisions [48|, we determine the averaged value of 
impact-parameter b = 3.15 fm for Woods-Saxon nuclear 
distribution and b = 2.2 fm for the hard-sphere distribu- 
tion. 



B. (l+3)d Ideal Hydrodynamical Expansion 

To take into account of both longitudinal and trans- 
verse expansion in the ideal hydrodynamic description of 
the bulk matter evolution, we use the output of calcu- 
lations by Hirano et at. [38l . f39| for Au+Au collision at 
the RHIC energy with Glauber collisions geometry and 
the average impact parameter b = 3.2 fm, corresponding 
to 10% most central events. The initial condition for the 
(l+3)d ideal hydrodynamic equations is fixed such that 
the final bulk hadro n sp ectra from the RHIC experiments 
are reproduced [H, |39| . We use the code provided by Hi- 
rano [49( that interpolates the numerical data from the 
hydrodynamic solution on a space-time grid for energy 
density, temperature, flow velocity and fraction of QGP 
phase. 

For our purpose, we have to infer parton and hadron 
number density from the energy density or temperature 
provided by the hydrodynamic solution, which uses a Bag 
model for the equation of state (EOS) [Elj. In this Bag 
model, the parton number and energy density in the QGP 
phase are, 



16^ + 36^7-^3, 



7T- 
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e QGP = (16^ + 36^JT 4 + B = 6iT* + B,(29) 

for three quark flavors. The Bag constant B = 
(247 Mev) 4 gives rise to a first-order phase transition 
from the QGP to a hadron resonance gas at the criti- 
cal temperature T c = 170 McV. 

In our calculation of medium modified parton fragmen- 
tation functions we consider jet-medium interaction in 
both partonic and hadronic phase throughout the evolu- 
tion of the bulk matter. Neglecting hadron correlation in 
the medium, the jet transport parameter will be propor- 
tional to hadron density and the soft gluon distribution 
within each hadron. There have been several phenomeno- 
logical studies [5lT[53j of jet quenching in deeply inelas- 
tic scattering off large cold nucleus as measured by the 
HERMES [54[ experiment. Within the same high-twist 
approach, the jet transport parameter is found [27] to be 
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FIG. 3: The temperature dependence of q/T 3 



qjv ~ 0.02 GeV 2 /fm at the center of the cold nucleonic 
matter in a large nucleus. In the (l+3)d ideal hydrody- 
namic model, the hadronic phase of the medium evolu- 
tion is considered as a hadron resonance gas, in which 
the jet transport parameter can be approximate as, 



qh = 



PN 



J2pm(t) + J2pb(t) 



M 



(30) 



where pn = Ho ~ 0.17 fm~ 3 is the nucleon density in 
the center of a large nucleus and the factor 2/3 accounts 
for the ratio of constituent quark numbers in mesons and 
baryons. The hadron density at a given temperature T 
and zero chemical potential is 
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\ T J 



E 



a (n^),(31) 



where rjh = ± for meson (M)/baryon (B). In the paper, 
we will include all hadron resonances with mass below 1 
GeV. Other hadron resonance gas models normally in- 
clude hadrons with mass up to 2 GeV which will only 
make some differences to our assumed qh close to the 
phase transition temperature. 

With the above approximation for hadron density and 
jet transport parameter in the hadron phase of bulk 
medium evolution, the jet transport parameter through- 
out the evolution of the medium can be expressed as, 



q{r,r) = q - 



JQGP 



p 



QGP 



(r ,0) 



(l-f) + q h (T,r)f. 



(32) 



where /(r, r) is the fraction of the hadronic phase at any 
given time and local position and qo denotes the jet trans- 
port parameter at the center of the bulk medium in the 
QGP phase at the initial time r . To illustrate the be- 
havior of the jet transport parameter in dense medium, 
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wc plot q/T 3 in Fig.[3]as a function of the temperature T 
according to the above assumptions, where we have used 
the initial value of go = 0.9^o'o4 GeV 2 /fm, as extracted 
by phenomenological study of experimental data later in 
this paper, at an initial temperature Tq = 369 MeV in 
the most central Au + Au collisions at RHIC. Appar- 
ently, we have neglected any additional temperature de- 
pendence beyond that in Eq. (|32[) . Pcrturbative calcula- 
tions [HI, HU in finite temperature QCD with resummed 
hard-thermal loops give rise to a logarithmic tempera- 
ture dependence of q/T 3 for a fixed value of the strong 
coupling constant a s . The same perturbative calcula- 
tions also lead to somewhat weak jet-energy dependence 
of the jet transport parameter which we will neglect in 
this phenomenological study. Therefore, one can consider 
the exacted jet transport parameter as an averaged value 
over the range of jet energies in the experimental data. 

In order to study the effect of phase transition in the 
(l+3)d hydrodynamic evolution of the bulk matter in 
jet quenching, we set / = in Eq. p"2")) and determine 
the parton density p = ai[(e — B)/b\} 3 / A according to 
Eq. |28|) throughout the bulk evolution, with e given by 
the (l+3)d hydrodynamic evolution at given space and 
time. Calculations with this prescription of space-time 
profile of the jet transport parameter will be labelled as 
no phase transition in this paper. We can also set q$ = 
to quantify the hadronic contribution to jet quenching. 

The jet transport parameter q in Eq. (|32[) is given in 
the frame where the medium is at rest. To take into 
account of the radial flow, one can simply make the fol- 
lowing substitute [H, H3] 



where p^ is the four momentum of the jet and m m is the 
four flow velocity in the collision frame, which is provided 
by the solution of the (l+3)d hydrodynamical equations 
[38l |39| . Such flow dependence effectively decreases (in- 
creases) the parton energy loss when the jet propagates 
along the same (opposite) direction as the flow. 

C. (l+3)d Parton Cascade 

To study the effect of non-equilibrium evolution of 
the bulk matter on jet quenching in this paper, we 
also consider a parton cascade model. The Boltzmann 
Approach of MultiParton Scatterings (BAMPS) model 
|44l 145^ solves the Boltzmann equation for on-shell glu- 
ons with pQCD based interactions, which include clastic 
scattering, brcmsstrahlung and its back reaction. The 
thermal equilibration process of gluons has been thor- 
oughly investigated within this parton transport model 
[43, [58[. BAMPS can in principle describe collective flow 
behavior j59Tj6l| and the nuclear modification factor Raa 
of jet quenching (62l | within a common framework. 

In this work we use the output from BAMPS [(30[ model 
for central Au+Au collisions with b — 3.4 fm at the RHIC 



energy. Initial gluons or minijets arc produced by binary 
nucleon-nucleon collisions with a Glauber geometry. The 
Lorentz contracted nuclei have a longitudinal extent of 
0.2 fm. Thus, most minijets are produced at t = 0.1 
fm/c when two nuclei fully overlap. The strong QCD 
coupling constant is set to be a constant, a s = 0.3. Glu- 
ons freeze out when the local energy density goes below a 
critical value, which is chosen as e c = 0.6 GeV/fm 3 . This 
corresponds to a critical temperature of T c = 175 MeV in 
a fully equilibrated gluon gas. In the present implemen- 
tation of BAMPS there are no hadronization and hadron 
cascade. Gluons propagate freely (free-streaming) after 
freeze-out. During this stage we regard one gluon as 
one hadron according to a parton-hadron duality picture. 
With these setups the final transverse energy and the el- 
liptic flow V2 from BAMPS calculations agree with the 
experimental data at RHIC (59[. 

We use Eqs. (|3"2"|) and (f3"3"|) to evaluate the jet transport 
parameter in the evolving bulk medium as described by 
the BAMPS model. The local gluon density and flow ve- 
locity are obtained by integrating the local gluon phase 
space distribution given by BAMPS. The hadron frac- 
tion / is either before gluon freeze-out or 1 after gluon 
freeze-out as a simplification for the hadronization pro- 
cess. By setting / = all the time, one effectively ne- 
glects the effect of phase transition and the difference 
in the jet transport parameter in QGP and hadronic 
phase. The hadron density will be estimated from a 
free-streaming hadron gas according to an assumption 
of hadron-parton duality. 



V. NUMERICAL RESULTS AND 
COMPARISONS 

With the three different models for the bulk matter 
evolution and the prescriptions for the evaluation of the 
jet transport parameter in the evolving medium as de- 
scribed in the last section, we can conduct a system- 
atic analysis of the experimental data on suppression 
of high pt hadron spectra within the higher-twist ap- 
proach to the medium modified jet fragmentation func- 
tions. Our aim in this work is to study the theoretical 
uncertainties associated with the dynamical evolution of 
the bulk medium and effects of jet hadron interaction 
in the hadronic medium on jet quenching. Since we as- 
sume that the jet transport parameter in the hadronic 
phase can be estimated from the value as determined in 
the cold nuclear matter in the DIS off nuclear targets 
[Eq. pop ], our specific problem in this case is to deter- 
mine the value of the jet transport parameter go at the 
center of the overlapped region in the central Au + Au 
collisions in the QGP phase at the initial time To [see 
Eqs. and l|32|)]. 

Shown in Fig. [Ha) are comparisons between the 
PHENIX experimental data on Raa for the neutral pion 
spectra in the 10% most central Au + Au collisions [§3[ 
at the highest RHIC energy yfs = 200 GeV/n and 
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FIG. 4: (color online) (a) Nuclear modification factor at midrapidity for the most 10% central Au+Au collisions at the 
maximum RHIC energy y/s = 200 GeV/n. The symbols are the PHENIX data taken from Ref. The curves are NLO 

pQCD calculations with three different values of goi~o in the (l+l)d Bjorken model with the Hard-Sphere (upper) and Woods- 
Saxon (lower) distribution function, (b) The corresponding \ 2 °f the fit as a function of <joTo- 



our NLO pQCD calculations with the medium modified 
fragmentation functions as given in the higher-twist ap- 
proach, for three different values of qo at the initial time 
tq = 0.6 fm/c, using (l+l)d Bjorken model of ideal hy- 
drodynamical evolution without a phase transition. The 
upper panel of Fig. [4] (a) shows the results with the 
hard-sphere nuclear geometry while the lower panel is 
for Woods-Saxon nuclear distribution. The difference 
between the calculated suppression factors with hard- 
sphere and Woods-Saxon nuclear distributions is about 
10 percent, due to the small difference in the thickness 
and overlap functions with these two nuclear distribu- 
tions. 

In the remainder of this paper, we will discuss the 
determination of initial jet transport parameter and its 
uncertainties associated with the modeling of the bulk 
medium evolution. We will carry out systematical \ 2 fits 
to the experimental data on the suppression factor Raa 
of neutral pion spectra from the PHENIX experiment 
[63| with our NLO pQCD calculations that include jet 
quenching in an evolving bulk medium. We will follow 
Ref. [HI to carry out a modified \ 2 analysis and ob- 
tain the best fit parameters go T o for each model of bulk 
medium evolution by minimizing 

x2 = ^ [Vi + e b a bi +w<r e - yff + ^ + ^ (JJ4) 
i=i Ui 



of the fits with respect to two parameters and e c , where 
t/i is the central value of experimental data with three 



types of uncertainties, <t,; (statistical plus uncorrelated 
systematics) , aj, t (correlated systematics) and a c (nor- 
malization uncertainties) for pr larger than 6.5 GeV/c, 
and y\ h is the calculated results at each pr point. Shown 
in Fig. @] (b) as two examples are the % 2 distributions 
as functions of go r o (to=0.6 fm/c) from the fits with 
the NLO pQCD calculation with (l+l)d Bjorken model 
for the bulk medium evolution. The best fit value is 
q T Q = 0.36GeV 2 and 0.39 GeV 2 for Woods-Saxon and 
hard-sphere nuclear distribution function, respectively. 
In the following calculations, we will always use the 
Woods-Saxon nuclear distribution function. 

Listed in Tabic Q] arc the best fit values of go To and the 
corresponding x 2 for jet quenching in various models for 
the bulk medium evolution. We will discuss and compare 
them in detail in the remainder of this section. 



A. (l + l)d Bjorken versus (l+3)d Ideal 
Hydrodynamical Model 

We first compare the results from (1+1 )d Bjorken and 
(l+3)d ideal hydrodynamic model with different options 
and study the effects of the transverse expansion, radial 
flow and contributions from jet quenching in the hadronic 
phase. Comparing the (l+l)d Bjorken and (l+3)d ideal 
hydrodynamical evolution without radial flow and phase 
transition it is clear that the transverse expansion cools 
the system considerably faster in the later stage of the 
evolution. This faster cooling due to the transverse ex- 
pansion leads to a faster decrease of the jet transport 
parameter in the (l+3)d hydrodynamic model than in 
the (l+l)d Bjorken model as shown in Fig.[S]by the dot- 
dashed [(l+3)d hydrodynamics] and solid lines [(l+l)d 
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Model 


flow 


PT q r (GeV 2 ) 


x 2 


(l+l)d Bjorken [W-S] 






0.36 


78 


(l+3)d Hydro 






0.42 


116 








0.48 


80 




V 




0.54 


106 


(<Zo = 0) 


V 


V 


0.78 


71 


BAMPS 






0.72 


95 








0.84 


84 




V 




0.96 


96 


(# = o) 


V 


V 


1.17 


86 



TABLE I: Initial jet transport parameter qoTo and the cor- 
responding x 2 - The check sign (-y/) denotes the inclusion of 
flow or phase transition (PT) in the calculations, respectively. 
The initial to is 0.6 fm/c for hydrodynamic and 0.3 fm/c for 
BAMPS model. 



Bjorken hydrodynamics]. To compensate for this faster 
cooling, one has to increase the value of the initial jet 
transport parameter q by 17% in the case of (l+3)d 
hydrodynamic evolution relative to the case of (l+l)d 
Bjorken expansion in order to fit the measured sup- 
pression factor Raa(pt) for neutral pions in the central 
Au + Au collisions at RHIC. 

Inclusion of the radial flow in the calculation of the jet 
transport parameter as in Eq. (|33[) will also change the 
effective jet quenching throughout the evolution of the 
bulk medium. Shown in Fig. [6] are the time evolution of 
the radial flow velocity in (l+3)d hydrodynamics (dashed 
line). The radial flow velocity in the center of the dense 
medium (left panel), where jet quenching is the strongest, 
remains small, at about 0.2, throughout hydrodynamic 
expansion. In the same figure, we also plot the fraction 
of the hadron phase (solid lines) as a function of time in 
the (l+3)d hydrodynamics. The period between / = 
and / = 1 indicate the duration of the phase transition 
or the mixed phase. We notice that at the beginning of 
the mixed phase the radial flow velocity starts to satu- 
rate and even decrease a little. This is mainly due to 
the mass effect. The energy, approximately ttitv, where 
niT is the transverse mass, should be the same during 
the phase transition. Therefore, the radial flow veloc- 
ity will decrease as particles become massive during the 
phase transition. Such decrease will be partially compen- 
sated by the continued collective expansion which tend 
to increase the flow velocity. As the phase transition is 
completed, the radial flow velocity should increase again 
by further transverse expansion in the hadronic phase 
and the resonances decays. 

The radial flow will reduce the effective jet transport 
parameter q if jets propagate along with it while increase 
the effective value of q if jets propagate against it. After 
averaging over the direction and initial position of jet pro- 
duction, the effective jet transport parameter is merely 
reduced by less than 14% [56|. As shown in Table HI this 



leads to an increase of go about 14% when one includes 
the effect of radial flow in the (l+3)d hydrodynamic evo- 
lution in fitting experimental data on Raa{pt)- 

In both (l+l)d Bjorken and two options of the (l+3)d 
hydrodynamic expansion as discussed above, we have ne- 
glected the difference of the jet transport parameters in 
QGP and hadronic phases. These options are labelled as 
no phase-transition (PT) by setting / = in Eq. (|32|) and 
using pQ GP = ai[(e - B)/^} 3 ^ in Eq. ([H throughout 
the evolution even in the mixed and hadronic phase. 

In the realistic scenario, jet transport parameter qh in 
the hadronic phase should be different from that in the 
QGP phase and we assume that qh is proportional to 
<7tv in the cold nuclear medium and the hadron number 
density as in Eq. (|30p . With this assumption, there is a 
discontinuity in q as a function of time (or temperature 
as in Fig. [3]) at the end of the mixed phase, as show by the 
dotted lines in Fig. [5l and the jet transport parameter in 
the hadronic phase in the later stage of the bulk medium 
evolution will be smaller than the scenario of no PT in 
the (l+3)d hydrodynamics. Therefore, one needs a larger 
value of the initial jet transport parameter q$, about 13%, 
to account for the measured suppression factor Raa(jPt) 
as compared to an evolving bulk medium with the same 
jet transport parameter throughout difference phases. 

To determine the contribution of jet quenching in the 
hadronic phase, we set q^ = in Eq. (|3"2"]l. assuming par- 
ton energy loss only happens in the QGP phase. In this 
case, the extracted qo will be about 44% larger (see Ta- 
ble []}. This implies that the hadronic phase contributes 
to about 44% of the total suppression of the high px 
hadron spectra. Such a large contribution is due to the 
lifetime of the mixed and hadronic phase until the ki- 
netic freeze-out which is longer than the lifetime of the 
QGP phase, though the hadron density is much smaller 
than the initial parton density in the QGP phase. A bet- 
ter understanding of the hadronic contribution to the jet 
quenching is, therefore, important for an accurate extrac- 
tion of the jet transport parameter in the initial phase of 
strongly interacting QGP. 



B. (l+3)d Parton Cascade versus (l+3)d Ideal 
Hydrodynamical Model 

In a parton cascade model such as BAMPS, the system 
takes times to reach thermal or partial thermal equilib- 
rium after the initial parton production time, tq = 0.3 
fm/c. During this period of parton equilibration, there 
is entropy production accompanied by rapid expansion. 
The net effect is a much faster decreasing of the effective 
temperature 0, E^, [65[ and the parton density during 
the early stage of the bulk medium evolution. The cor- 
responding jet transport parameter in the equilibrating 
gluonic matter in BAMPS calculation therefore decreases 
much faster than that in both (l+l)d and (l+3)d ideal 
hydrodynamic models in the early stage of evolution as 
shown by the dashed lines in Fig. [5] It is smaller by 
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FIG. 5: (color online). Time evolution of the scaled jet transport parameter at position of radius of 2 fm (left panel) and 4 fm 
(right panel) for various models of the bulk matter evolution. The flow effect in Eg . (|33[l is not included. 
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FIG. 6: (color online). The radial flow velocity in the (l+3)d hydrodynamical model and the parton cascade at position of 
radius of 2 fm (left panel) and 4 fm (right panel). The solid curves show the fraction of the hadron phase in the hydrodynamical 
model. 



more than a factor of 2 before the hadronization. This 
leads to about 77% increase in the initial value of qq 
that is needed to account for the measured suppression 
of hadron spectra (see Tabic |T|). Note that in our calcula- 
tions jet quenching starts at the thermalization tq = 0.6 
fm/c in the ideal hydrodynamic model. However, it starts 
at To = 0.3 fm/c in the BAMPS model immediately af- 
ter the initial gluons are produced and before thermal or 
partial thermal equilibrium is reached. Such jet medium 
interaction and parton energy loss during the thermaliza- 



tion stage in the BAMPS model will contribute to some 
additional, though a small fraction, suppression of the 
final hadron spectra. 

In BAMPS, the hadronization of the gluonic matter 
happens when the local gluon energy density is below 
e c = O.GGcV/ fm 3 and parton-hadron duality is used for 
the phase-space density of hadrons, which are assumed 
to freeze out immediately. Such a sudden hadronization 
and freeze-out scenario implies the change of the hadron 
fraction / from to 1 at this critical density in the eval- 
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uation of the jet transport parameter in Eq. (|32[) . The 
lack of the mixed phase in BAMPS, which lasts for about 
5 — 6 fm/c in the (l+3)d ideal hydrodynamic evolution 
as indicated by the duration for the hadron fraction to 
reach / = 1 in Fig. [6j leads to smaller value of the jet 
transport parameter after the hadronization. Comparing 
the extracted values of go from the phenomenological fit 
with and without (q^ = 0) hadronic phase as shown in 
Table U the hadronic interaction in BAMPS model con- 
tributes to about 22% of the jet quenching effect, which 
is about a factor 2 smaller than that in the (l+3)d ideal 
hydrodynamical model. 

The large viscous effect in the parton cascade model 
generally leads to larger radial flow velocity than that 
generated by an ideal hydrodynamic expansion as shown 
in Fig. |6l The free-streaming of hadrons immediately 
after the hadronization also leads to the further increase 
of the flow velocity. Such larger radial flow velocity in 
the BAMPS model reduces the effective jet quenching 
parameter and therefore leads an increase of about 16% 
in the extracted initial qo . 

The decrease of the radial flow velocity at late times 
seen in Fig. [5] is a numerical artifact, which comes from 
the chosen geometry of spatial elements, where particles 
are collected to calculate the radial flow velocity. Longi- 
tudinal length of each spatial element corresponds to a 
space time rapidity window \rj\ < 0.25, while the trans- 
verse length is fixed to be 0.3 fm. At late times the 
longitudinal length is larger than the transverse one, and 
particles with smaller transverse velocity (larger longi- 
tudinal velocity) make up the main part of remaining 
particles in the spatial elements. This is the reason for 
the observed decrease of the radial flow velocity. If we 
choose the longitudinal length being equal as the trans- 
verse one at the late times, we would see a continuous 
increase of the flow velocity due to the free-streaming, 
before all particles go out of the spatial element. This 
is, however, a negligible effect on the extraction of qo, 
because the hadron density is tiny at late times. 



VI. CONCLUSIONS 

Using medium-modified fragmentation functions from 
the high-twist approach to multiple parton scattering we 
have calculated single inclusive hadron spectra in NLO 
perturbative QCD parton model. The medium modifi- 
cation to the fragmentation function is very similar in 
form to the DGLAP correction due to vacuum gluon 
brcmsstrahlung, except that the medium modified split- 
ting function contains information about the properties 
of the medium through a path-integrated jet transport 
parameter q. Therefore, the medium modified fragmen- 
tation function will explicitly depend on the space-time 
evolution of the bulk medium. We have focused our at- 
tention on the dependence of the medium modified frag- 
mentation function on the space-time evolution of the 
bulk medium and their influence on the suppression of 



the single hadron spectra at large pr- 

We specifically considered the effects of transverse ex- 
pansion, radial flow, jet quenching in hadronic phase and 
non-equilibrium effect with three different models for the 
bulk matter evolution in central Au + Au collisions at 
the RHIC energy: (l+l)d Bjorken hydrodynamic model, 
(l+3)d ideal hydrodynamical model by Hirano [H, |39| 
and (l+3)d BAMPS HHH] parton cascade model. Since 
the jet transport parameter is proportional to the parti- 
cle gluon distribution density of the medium, we have 
assumed that it will be proportional to the local particle 
density in both QGP and hadronic phase during the bulk 
matter evolution. The coefficient of the particle density 
dependence, which is related to soft gluon distribution 
per particle, will be different in the QGP and hadronic 
phase. Assuming a similar density dependence of the jet 
transport parameter in the hadronic phase as that in a 
large cold nuclei, but rescaled by the local hadron density 
and the number of constituent quarks, the only free pa- 
rameter in the space-time profile is the initial value of the 
jet transport parameter qo at initial time To at the cen- 
ter of the QGP medium in central heavy-ion collisions. 
We have extracted values of go a t t m each model of 
bulk evolution with different options by fitting the ex- 
perimental data on the suppression factor Raa{pt) for 
high pt neutral pions and studied the uncertainties in 
the extracted value of go- 

By comparing the values of go extracted with different 
options in models for the bulk evolution, we also found 
that the transverse expansion leads to fast cooling of the 
hot medium and therefore reduces jet quenching in the 
later stage of the evolution by about 17%. Radial flow 
also reduces the effective jet transport parameter by 15%. 
The phase transition from QGP to hadronic medium fur- 
ther reduces the effective jet transport parameter in the 
hadronic phase and leads to about 13% reduction of jet 
quenching effect. Within our model for the form of jet 
transport parameter due to jet-hadron interaction, the 
overall contribution to jet quenching from hadronic phase 
is about 22-44%, depending on the model for evolution 
of the hadronic phase. 

By comparing (l+3)d ideal hydrodynamic and the 
BAMPS parton cascade model, we found that the non- 
equilibrium evolution in the early stage of parton cascade 
leads to faster decrease of the jet transport parameter 
with time and therefore affects the overall jet quenching 
the most by 77%, part of which is caused by the lack of 
mixed phase in the parton cascade model. Such uncer- 
tainties can be further reduced by combining parton cas- 
cade and hadron cascade with a model of of hadronization 
that can reproduce the EOS during the phase transition 
as given by lattice QCD calculation. 
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